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We study the diversities in the properties of the neutron stars arising due to the different choices for 
the cross-coupling between various mesons which governs the density dependence of the nuclear sym¬ 
metry energy in the extended relativistic mean-field(RMF) model. For this purpose, we obtain two 
different families of the extended RMF model corresponding to different non-linear cross-coupling 
term in the isovector part of the effective Lagrangian density. The lowest order contributions for the 
S mesons are also included. The different models within the same family are so obtained that they 
yield wide variation in the value of neutron-skin thickness in the ^°®Pb nucleus. These models are 
employed to compute the neutron star properties such as, core-crust transition density, radius and 
red shift at canonical mass (I. 4 M 0 ), tidal polarizability parameter, and threshold mass required for 
the enhanced cooling through direct Urea process. Most of the neutron star properties considered 
are significantly different (10%-40%) for the different families of models at a smaller neutron-skin 
thickness (~ 0.15 fm) in the ^°®Pb nucleus. 
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I. INTRODUCTION 

The neutron stars are believed to be composed of 
highly asymmetric matter, predominantly, the neutrons. 
Of course, a small admixture of protons, electrons and 
muons are also present to maintain the /3-equilibrium and 
charge neutrality. The various theoretical models have 
conjectured the possibility of existence of exotica, like, 
hyperons, Bose condensates and quarks in the core of the 
neutron stars [in. The precise knowledge of the masses 
and the radii of the neutron stars are crucial in deter¬ 
mining their compositions. Recently, significant progress 
has been made along this direction [^. The masses of 
PSR J1614-2230 0 and PSR J0348-b0432 [nj are mea¬ 
sured to be 1.97 ± O.O4M0 and 2.01 ± O.O4M0, respec¬ 
tively. These measurements impose the lower bounds on 
the maximum mass of the neutron stars that a theoretical 
model must yield. In the absence of the exotic degrees of 
freedoms, the criterion of maximum mass to be ^ 2Mq 
is readily satisfied by the theoretical models. The ob¬ 
servational constraints on the maximum mass of neutron 
stars do not completely rule out the possibility of exis¬ 
tence of the exotic degrees of freedom, but, the threshold 
transition densities for their appearance are pushed to 
2.5-3.5 times the nuclear saturation density |12 h 14|| . The 
radii of neutron stars are known only poorly. The values 
of radii are quite sensitive to the assumed composition 
of the atmosphere. The radius Ria = 10.7 — 13.1km, 
for the neutron star with the canonical mass of 1.4M0, 
is found to be consistent with the observational analysis 
and the host of experimental data for finite nuclei [15l | . 
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The equations of state (EOSs) for the /3-equilibrated 
baryonic matter, employed to determine the bulk prop¬ 
erties of neutron stars, are usually constructed using the 
energy density functionals derived from the Skyrme type 
effective forces or from an effective Lagrangian density 
associated with the relativistic mean field (RMF) model. 
Often, the energy density functionals are optimized us¬ 
ing some selected experimental data on a key properties 
of the finite nuclei. Occasionally, pseudo-data on nuclear 
and neutron matter are also used in the optimization pro¬ 
tocols [mil. The EOS for the nuclear matter at a given 
density and the asymmetry can be viewed for simplicity 
as, 

e{p,I) = e{p,0) + S{p)I^ (1) 

where, e{p,I) is the energy per nucleon, p = Pn + Pp, 
I = (Pn — Pp)lP is the asymmetry with and pp being 
densities for the neutrons and the protons. The e{p, 0) 
is the EOS for the symmetric nuclear matter and S{p) is 
the symmetry energy coefficient at a density p. The EOS 
of the symmetric nuclear matter for the densities up to 
4.5 times saturation density (po = 0-16 fm~^) is derived 
within the reasonable limits by combining the experi¬ 
mental data for the finite nuclei with the collective flow 
and kaon production data in heavy-ion collisions (Tol - f^ . 
The poorly known density dependence of S{p) is the ma¬ 
jor source for uncertainty in the EOS for the asymmetric 
matter . The value of S{p) is reasonably constrained 
only around the saturation density by the bulk proper¬ 
ties of the finite nuclei. The understanding of density 
dependence of the S{p) is crucial as it controls the radii 
of neutron stars, the thicknesses of their crusts, the rate 
of cooling of neutron stars, and the properties of nuclei 
involved in r-process nucleosynthesis. The density depen¬ 
dence of the symmetry energy around the saturation den- 
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sity appears to be well correlated with the neutron-skin 
thickness Ar^p in a heavy nucleus which can be experi¬ 
mentally measured. The Arnp is the difference between 
the rms radii for density distributions of the neutrons 
and protons in a nucleus. Recently [13, the correla¬ 
tions of the neutron-skin thickness in ^°®Pb nucleus with 
several bulk properties of neutron stars have been exam¬ 
ined for the TOV-min and FSU type models. The energy 
density functional for the TOV-min corresponds to the 
Skyrme type effective force and that for FSU is based 
on the extended RMF model. The correlation between 
neutron-skin thickness in ^°®Pb nucleus and the neutron 
star radius Ria for TOV-min is noticeably smaller than 
the one obtained for the FSU model. Consequently, for 
the case of TOV-min the properties of neutron stars can 
have larger variations at a fixed neutron-skin in ^°®Pb 
nucleus. This result is in concordance with the large un¬ 
certainties in the high density behaviour of the symmetry 
energy for the Skyrme type energy density functionals 


24|. 


The tighter correlations of neutron-skin thickness in 
208pb nucleus with the several properties of the neutron 
stars within the RMF models seem to be stemming from 
the lack of freedom in the isovector part of the effective 
Lagrangian density associated with these models. Most 
of the RMF models describe the density dependence of 
the symmetry energy either in terms of the coup ling of 
the isovector-vector p mesons with the nucleons |25l | or 
by including only an additional cross-coupling of isosclar- 
scalar a or isosclar-vector w with the isovector-vector p 
mesons [H, [23| . The a — p oi oj — p cross-coupling al¬ 
lows one to vary the neutron-skin thickness over a wide 
range without significantly affecting the quality of fit to 
the bulk properties of the nuclei such as the total bind¬ 
ing energy and the charge radii 23] . The differences in 
the high density behaviour of the symmetry energy and 
their consequences on the properties of the neutron stars 
arising due to the use of different cross-couplings have 
never been studied in detail. Further, the inclusion of 
isovector-scalar S mesons can modify the behaviour of 
the symmetry energy at high densities [23 - l3lj |. In prin¬ 
cipal, the contributions from the various cross-couplings 
and the 5 mesons should be included in a single RMF 
model. It has not been done so far due to lack of accu¬ 
rate experimental data on the finite nuclei and the neu¬ 
tron stars which govern the isovector part of the RMF 
model. A comprehensive study of variations in the prop¬ 
erties of the neutron stars as a function of the neutron- 
skin thickness within the RMF models corresponding to 
different choices for the cross-coupling terms with and 
without the inclusion of the S mesons may be highly de¬ 
sirable. Such investigation would enable one to under¬ 
stand to what extent the inclusion of contributions from 
these various cross-couplings and the S mesons within a 
single RMF model are necessary to describe simultane¬ 
ously the neutron-skin thickness in ^°®Pb and the various 
neutron star properties. 

In the present work, we would like to study the diver¬ 


sities in the properties of the neutron stars arising purely 
due to the uncertainties in the isovector part of the ef¬ 
fective Lagrangian density which governs the density de¬ 
pendence of the nuclear symmetry energy in the extended 
RMF model. Towards this purpose, two different fami¬ 
lies of extended RMF models are obtained which mainly 
differ from each other in the choice for the cross-coupling 
term in the isovector part of the effective Lagrangian den¬ 
sity. One of the families of models includes a — p cross¬ 
coupling while the other includes uj — p cross-coupling 
term in addition to the various linear and non-linear in¬ 
teraction terms already present in the commonly used 
RMF models. The contributions due the coupling of the 
S mesons to the nucleons are also considered. The var¬ 
ious coupling constants are so varied that they produce 
wide variations in the neutron-skin thickness in ^'^^Pb nu¬ 
cleus without affecting significantly the binding energies 
and charge radii of the finite nuclei. The various neu¬ 
tron star properties considered are the core-crust transi¬ 
tion density, radius for the neutron stars with canonical 
mass, the tidal polarizability parameter and the thresh¬ 
old mass required for the enhanced cooling through direct 
Urea process. Some of these neutron stars properties at 
a fixed neutron-skin thickness differ significantly for two 
different families of the models. 

The paper is organized as follows. In Sec. II, we briefly 
outline the form of the effective Lagrangian density for 
the extended RMF model. The procedure adopted to ob¬ 
tain different parameterizations for two different families 
of the model are described in Sec. III. We also present in 
this section the results depicting the relationship between 
different coupling constants which govern the isovector 
part of the effective Lagrangian density. In Sec. IV, the 
various properties of the neutron stars obtained for the 
different families of models are compared for fixed values 
of neutron-skin thickness in ^°®Pb nucleus. The main 
conclusions are presented in Sec. V. 


II. THEORETICAL FRAME WORK 

In the RMF models, nucleons interact through the ex¬ 
change of isoscalar scalar tr, isoscaler vector w and isovec¬ 
tor vector p mesons. The effective Lagrangian density for 
the RMF model usually includes the cubic and quartic 
order non-linear self-interaction terms for the a mesons 
in addition to the linear terms for the cr, w and p mesons 
which describe their interactions with the nucleons. The 
non-linear self-interaction terms for the a mesons are 
added to yield reasonable values for the empirical prop¬ 
erties of symmetric nuclear matter. Further, the RMF 
models are extended by including various cross-couplings 
terms for these mesons and self-interaction terms for the 
Lo and p mesons. The a — p and uj — p cross-coupling 
terms enables one to vary the density dependence of the 
symmetry energy coefficient and the neutron skin thick¬ 
ness in heavy nuclei over a wide range without affecting 
the other properties of finite nuclei [H, js^j • Most of the 
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RMF models do not include the contributions from the 
isovector-scalar <5 mesons. The bulk properties of the fi¬ 
nite nuclei like binding energies and radii are not very 
sensitive to the presence of the S mesons. However, ear¬ 
lier investigations have stressed the need to include the 
contributions from the <5 mesons for proper descri ptio n 
of the highly asymmetric dense matter. l34l - [^ . 

The effective Lagrangian density which includes the 
lowest order contribution from the 6 mesons together 
with the various non-linear cross-coupling and self¬ 
interaction contributions already present in the extended 
RMF model, can be written as, IstI-I^. 


£, — M “b + c,p + + c, aujp ; (2) 

where the Lagrangian Cj\fM describing the interactions 
of the nucleons through the mesons is. 


= X! ^ - ga(J - gsSr) 

J=n,p 

(3) 

Here, the sum is taken over the neutrons and protons 
and T are the isospin matrices. The Lagrangian density 
Ci for i = a,u}, p and S can be written as. 


^<7 = \idpCrd^cr-mlaA--^9<ymla^-^A^glrnla^, 

( 4 ) 








( 6 ) 


Cs = \{dp6d>^6-mjSA. (7) 

The u!^'' , p^'^ are field tensors corresponding to the oj and 
p mesons, and can be defined as — d'^uj^ and 

ppv _ QPpV _ Qv The cross interactions of cr, w, and p 
mesons are described by which can be written as. 


= ^9crmZ<7UJpUjf^ + -^^g^mZa^ujpUJ^ 
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m 


+ 


mp 


J^P^rUpUJpUJ'^PpP'- 


One also needs to include the contributions from the elec¬ 
tromagnetic interaction in the case of finite nuclei. The 
Lagrangian density Cem for the electromagnetic interac¬ 
tion can be written as. 


C 


em 


- e^pjpAp'i’p, 


(9) 


where, A is the photon filed and = d^A" — d''A^. 
The equation of motion for nucleons, mesons and photons 
can be derived from the Lagrangian density defined in 
Eq. (j^j). T'he contributions from Eq. m are included 
only for the case of finite nuclei. 

It is clear from dS]) that there are five cross-coupling 
terms; two of them are the cubic order terms correspond¬ 
ing to the (T — OJ and a — p cross-couplings and the re¬ 
maining are the quartic order terms. The contribution 
from the a — oj cross-couplings and self coupling of oj 
mesons play important role in varying the high density 
behaviour of the EOSs and also prevents instabilities in 
them SEaEIi. The contributions of the self-coupling 
of p mesons are not considered in Eq. m, since, expec¬ 
tation value of the p meson field is order of magnitude 
smaller than that for the oj meson field . The inclusion 
of the p meson self interaction can affect the properties of 
the finite nuclei and neutron stars only very marginally 
(4ll | . Of the particular interest in the present work are the 
cross-coupling terms involving p meson field which con¬ 
tributes to the isovector part of the effective Lagrangian 
density in addition to the usual linear couplings of the 
p and 5 mesons to the nucleons. We shall mainly fo¬ 
cus on the lowest order a — p and oj — p cross-couplings 
whose strengths are determined by the values of Pp and 
P 2 p- The quartic order a — p cross-coupling strength pip 
is set to zero. The values of pp or p 2 p can be appro¬ 
priately adjusted to yield wide variations in the density 
dependence of the symmetry energy coefficient and the 
neutron skin thickness in heavy nuclei without affecting 
the other properties of finite nuclei [H, HI, . 


III. MODEL PARAMETERS 

Two different families of the extended RMF models, 
named hereafter as Fp and Ejp, are obtained. These 
families differ from each other in the choice for the cross¬ 
coupling term in the isovector part of the Lagrangian 
density. The isovector part of the Lagrangian density for 
the Fp{F 2 p) family is governed by the coupling parame¬ 
ters gp, gs and Pp{p 2 p)- The parameters gp and gs denote 
the strengths for the coupling of the p and 6 mesons to 
the nucleons, respectively. The parameter pp and p 2 p de¬ 
note the strength of the a — p and oj — p cross-couplings as 
can be seen from Eq. [51 The remaining parameters which 
correspond to the isoscalar part of the Lagrangian den¬ 
sity and the mass of the a, oj and p mesons are kept fixed 
to that of the BKA22 model [i^. The BKA22 model 
has been identified to satisfy various constraints related 
to symmetric nuclear matter, pure neutron matter, sym¬ 
metry energy, and its derivatives [d^. 

The different parameterizations of Fp{F 2 p) families are 
obtained by varying appropriately the values of gp, gs and 
Vp{V 2 p)- For a given value of gs and Pp{p 2 p), the value of 
gp is always adjusted to yield appropriate binding energy 
for the ^°®Pb nucleus. Once the values of gs, gp and pp or 
P 2 p are known, the properties of the nuclear matter and 
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TABLE I: The parameters of the isovector part of the La- 
grangian density for some representative sets for the Fp and 
F2, families of models. In the bottom part, the values for 
the symmetry energy coefficient at the saturation density J, 
symmetry energy slope parameter L, effective mass for the 
protons and neutrons and their differences are also presented. 
All these quantities are in MeV. 


Parameter 

J 

P 

F2p 

SETl 

SET2 

SETS 

SET4 

»7p 

4.0 

4.0 

0.0 

0.0 

V2p 

0.0 

0.0 

17.5 

17.5 

gs 

0.0 

8.0 

0.0 

8.0 

9p 

13.033 

21.863 

11.051 

18.556 

J 

33.3 

30.9 

33.0 

30.2 

L 

79.0 

62.5 

65.0 

23.3 

m; 

570.1 

630.9 

577.8 

652.4 

m: 

570.1 

495.4 

577.8 

515.8 

AMp*„ 

0.0 

135.5 

0.0 

136.6 


TABLE 11: The values of the total binding energy {E ) in 
MeV, charge radii (rc), neutron radii (rn) and neutron-skin 
thickness Arnp in fm for a few asymmetric spherical nuclei 
obtained for SETl - SET4 parameters. 


Nucleus Property 

Fp 


F2p 

SETl 

SET2 

SETS 

SETl 

'‘“Ca 

E 

-415.75 

-415.40 

-415.81 

-415.61 


Tc 

3.468 

3.484 

3.465 

3.477 


rn 

3.575 

3.544 

3.574 

3.535 



0.201 

0.153 

0.202 

0.151 

'^"Sn 

E 

-1102.49 - 

1100.37 

-1102.69 

-1100.99 


Tc 

4.736 

4.759 

4.727 

4.731 


Tn 

4.952 

4.901 

4.934 

4.873 



0.284 

0.210 

0.286 

0.210 


E 

-1637.07 - 

1637.08 

-1637.06 

-1637.05 


Vc 

5.545 

5.566 

5.535 

5.538 


Vn 

5.706 

5.659 

5.699 

5.631 



0.219 

0.151 

0.221 

0.151 


the finite nuclei can be computed. We vary the values of 
gs over a wide range from 0 to 8. The values of r]p and r) 2 p 
are varied in the range of 0 — 12 and 0 — 60, respectively. 
For r/p > 12, the stable solutions of the field equations for 
the mesons could not be obtained. We have constructed 
22 different parameterizations of the Fp and 41 different 
parameterizations of the F 2 p families. 

The Fig. [1] displays the relationship between the var¬ 
ious parameters of the isovector channel for the Fp and 
F 2 p families of the models. It can be readily seen that 
the value of gp, required to reproduce the binding energy 
for the ^*^®Pb nucleus, increases with gs, rip and rj 2 p- In 
other words,the equation of state at least for the densi¬ 
ties relevant for the finite nuclei becomes softer with the 
increase in the gs, rjp and ri 2 p, which is compensated by 
increasing the value of gp to reproduce the binding en¬ 
ergy for the ^'’®Pb nucleus. In the Fig. [2l the relationship 
of the parameters gs and rip{ri 2 p) with the symmetry en¬ 
ergy coefficient at the saturation density (J = S{po)) for 
the Fp{F 2 p) is displayed in terms of the contour plots. 
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FIG. 1: (Color online) Colour coded contours in the gs — 
r;p (upper panel) and gs — 772 P (lower panel) planes correspond¬ 
ing to the Fp and F 2 p families, respectively. The value of gp 
are colour coded according to the scale one the right side. 


Similarly, the results for the Ar^p in the ^°®Pb nucleus 
are plotted in Fig. In general, the values of J and 
Aunp decreases with increasing gs, Pp or p 2 p- The F 2 p 
model yields larger variations in Ar^p. The large values 
for pp are not favored, as a result the Fp family can yield 
very small values of Aunp only with the inclusion of the 
S mesons. In Table [J the values of the parameters for 
four representative sets corresponding to the Fp and F 2 p 
models are listed. The SETl and SET2 belong to the 
Fp family, while, SETS and SET4 are for the F 2 p family. 
The SETl and SETS do not include the contributions 
from the S mesons {gs = 0). The SET2 and SET4 cor¬ 
respond to the highest value of the (5-nucleon coupling 
strength {gs = 8), otherwise, they are very much simi¬ 
lar to the SETl and SETS, respectively. These different 
sets are so chosen that the comparison of the properties 
of the neutron stars resulting from them would give us 
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FIG. 2: (Color online)Same as Fig. [T] but, for the symmetry 
energy at the saturation density (J = S{po)) fixed along the 
contour. The values of J are in MeV. 

a crude estimate about the effects of 6 meson as well as 
the a — p and uj — p cross-couplings. In the bottom part 
of Table H the values of the symmetry energy coefficient 
at the saturation density (J = S{po)), symmetry energy 
slope parameter {L = ), proton and neutron 

effective masses and their differences are also presented. 
The effective masses are obtained at the maximum asym¬ 
metry, i.e., the pure neutron matter. In Table m we 
present some bulk properties of a few asymmetric spher¬ 
ical nuclei. The various bulk properties for these nuclei 
are relatively better reproduced for the SETl and SETS 
parameters which corresponds to Aunp ~ 0.22fm in the 
208pb nucleus, since, this value of Aunp is almost the 
same as that of the base model BKA22. It may be noted 
that the Aunp = 0.22 fm for the SETl and SETS, but, 
they belong to different families. Similarly, SET2 and 
SET4 represent different families, with Ar^p = 0.15 fm. 

Let us now take a look at the density dependence of 


FIG. 3: (Color onlme)Same as Fig. [U but, for the neutron- 
skin thickness Arnp in the ^*’®Pb nucleus fixed along the con¬ 
tour. The values of Arnp are in fm. 


the symmetry energy for the different parameterizations 
of the Fp and F 2 p families corresponding to equal values 
of Aupp. In Fig. |4l we display the variations of symmetry 
energy as a function of density for different parameter¬ 
izations as indicated by SETl, SET2, SETS and SET4. 
The high density behaviour for the symmetry energy is 
stiffer for the Fp family as can be easily verified by com¬ 
paring the results for the SETl and SET2 with those for 
the SETS and SET4, respectively. Further, by comparing 
the results for the SETl with SET2 or those for SETS 
with SET4, it can be concluded that the inclusion of the 
S mesons softens the symmetry energy at low densities 
while makes it stiffer at higher densities. The results 
depicted in Fig. 2] provide evidences a priori about the 
possibilities of the differences in the properties of neu¬ 
tron stars at a fixed Ar^p across the different families 
of the models, due to the differences in the high density 
behaviour of the symmetry energy. 



























IV. NEUTRON-SKIN THICKNESS AND 
PROPERTIES OF NEUTRON STARS 
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We wish to study the differences in the properties of 
neutron stars for the Fp and F2p families of the models 
at fixed values for the neutron skin thickness Arnp in the 
208pb nucleus. In particular, attention is given to the 
study of such differences at Arnp = O.lSfm in nu¬ 

cleus. This value of Arnp is consistent with 0.156lo.o2i 
fm [ 4 ^ and 0.168±0.022 fm as extracted from the ex¬ 
perimental data on the dipole polarizability for ^°®Pb nu¬ 
cleus. A ver y re cent measurement of coherent pion photo¬ 
production also corresponds to Arnp = 0.15 ± 0.03 
fm in ^*^®Pb nucleus. However, these measurements do 
not conclusively yet rule out the larger values for Arnn, 
since, the Lead Radius Experiment (PREX) [H, IH, 
has recently measured Ar„p = 0.33iliQ j^g fm in ^°®Pb nu¬ 
cleus via parity-violating electron scattering which pro¬ 
vides the first purely electroweak, almost model indepen¬ 
dent estimate. Our purpose is to give quantitative esti¬ 
mates about the extent to which the various properties 
of the neutron stars might vary across the different fam¬ 
ilies of the models for a plausible value of neutron-skin 
thickness. The various properties for the neutron stars 
considered are the core-crust transition density, radius, 
red-shift, the threshold mass required for the enhance 
cooling through the direct Urea process and the tidal po¬ 
larizability parameter. The comparison of results for the 
Fp and F^p families of the models would enable us to un¬ 
derstand the role of different cross-coupling terms. We 
shall also assess the effects of the 5 mesons by comparing 
the results obtained with and without its inclusion in the 
same family of the models. 

The values of the neutron-skin thickness in a heavy nu¬ 
cleus, according to the Droplet Model [H^, are strongly 
correlated with the symmetry energy slope parameter L. 
The dependence of L on Ar„p in ^°®Pb for the Fp and 
F 2 p families of the models are displayed in Fig. [S] It may 
be pointed out that the similar values of Arnp can be ob¬ 
tained within a given family by varying appropriately the 
values of coupling parameters gs and r]p or r] 2 p(see also 
Fig. [3]). The solid and the hollow symbols represent 
the results obtained with and without the contributions 
from the <5 mesons, respectively. The values of Arnp are 
well correlated with L within a given family of the mod¬ 
els irrespective of the contributions from the S mesons. 
However, the values of L for the two families of the mod¬ 
els differ significantly at smaller Arnp. This difference 
gradually disappears as Arnp increases. 

The values of core-crust transition density pt and the 
corresponding pressure Ft as a function of Arnp obtained 
for Fp and F 2 p families of models are plotted in Fig. 
E The values of pt are obtained using a method based 
on the relativistic random-phase approximation [5ll - l53 |. 
This method uses the fact that the uniform matter in its 
ground state at sufficiently low densities becomes unsta¬ 
ble to small density fluctuations. The values of pt are 
correlated with the Arnp within a same family irrespec- 



FIG. 4: (Color online) The density dependence of symmetry 
energy S{p) for some representative cases of Fp and F2p fam¬ 
ilies of the models. The labels SETl and SET2 correspond to 
the two different parameterizations for the Fp family, whereas, 
the SETS and SET4 correspond to the F 2 p family. The SETl 
and SETS are associated with Arnp = 0.22fm and the SET2 
and SET4 yield Arnp = 0.15fm (see also Tables III and [Till . 
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FIG. 5: (Color online) The variations of symmetry energy 
slope parameter L with neutron-skin thickness Arnp in the 
^°®Pb nucleus for the Fp and F2p families of models. The 
solid and hollow symbols represent the results obtained with 
[gs yf 0) and with {gs — 0), respectively. 


five of the contributions from the 6 mesons. But, this 
correlation seems to be some what model dependent— 
the values of pt for both the families of models at a fixed 
Arnp are not the same. In particular, the pt is signif¬ 
icantly larger for the F 2 p family at smaller Arnp. The 
transition pressure Pt is not very well correlated with 
the Arnp. Initially, the Pt increases with Arnp and it 
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FIG. 6: (Color online) The variation of core-crust transition 
density and the corresponding pressure with the neutron-skin 
thickness Arnp in nucleus for the Fp and Fbp families 

of the extended RMF models. 



FIG. 7: (Color online) The mass-radius relationship for the 
Fp and F2p families of the models. The solid squares represent 
the masses and the corresponding radii for the neutron stars 
with the central density to be Spo- 

decreases for higher values of Arnp. 

Once, the core-crust transition density is determined, 
the EOS for various density ranges as required for the 
computation of the properties of the neutron stars can 
be constructed. The EOS data for the density p ~ 
4.8 X 10“® — 2.6 X lO^^^fm”^ corresponding to the outer 
crust region are taken from Ref. [^. The EOS for the 
inner crust is obtained by assuming a polytropic form 
P(e) = a -I- where P and e are the pressure and 

energy density respectively. The constants a and b are 
determined in such a way that the EOS for the inner crust 
matches with that for the inner edge of the outer crust 


at one end and with the edge of the core at the other 
end. The EOSs for the core region, p > pt, are obtained 
within the RMF model by using the different parameter- 
izations of the Fp and F 2 p families. The core region is 
assumed to be composed of neutrons, protons, electrons 
and muons. The chemical potentials for various particle 
species at a given baryon density are obtained by impos¬ 
ing the /3-equilibrium and charge neutrality conditions. 
We use these EOSs to compute the properties of static 
neutron stars by integrating the Tolman-Oppenheimer- 
Volkoff (TOV) equations [^. In Fig. Qwe display the 
mass-radius relationship for the sequences of static neu¬ 
tron stars obtained for the Fp and P 2 p families of models. 
The solid and the dashed lines depict the results obtained 
with and without the inclusion of the contributions from 
the S mesons, respectively. The solid squares represent 
the masses and the corresponding radii for the neutron 
stars with the central density to be 3po. The different 
EOSs obtained for a given family of the models differ 
mainly in the high density behaviour of the symmetry 
energy. This leads to the variations in the mass-radius re¬ 
lationship for the neutron stars within the same family of 
the models.The maximum mass Mmax = 1-95 — 2.02Mq 
and 1.91 — 1.98Mq and the radii i?i ,4 = 13.3 — 15.4km 
and 12.3 — 14.9km for the Fp and F 2 p families, respec¬ 
tively. The value of Mmax for the Fp family is consistent 
with the recent mass measurements M = 1.97 ± 0.04Mq 
for PSR J1614-2230 [ 13 ] and M = 2.O1±O.O4M0 for PSR 
J0348-I-0432 [Tl|, but, the values for i?i ,4 is marginally 
away from the i?i ,4 = 10.7 — 13.1km as extracted in Ref. 
[l5| . For the F 2 p family, values for Mmax are barely con¬ 
sistent with the recent measurements, but, i?i ,4 is consis¬ 
tent with the ones extracted in Ref. [l3| . The inclusion of 
the 5 mesons yields higher values for the maximum mass 
for the neutron stars within a family. A more realistic 
estimation for the effects of 5 mesons on the maximum 
mass of the neutron stars requires the inclusion of vari¬ 
ous exotic degrees of freedom. Since, the density at the 
center of the neutron star with the maximum mass for 
our EOSs is significantly larger than Spo- 

We now compare the various properties of the neutron 
stars at fixed values of Arnp obtained for the Fp and 
F^p families. Before embarking on our discussion, it may 
be reminded that the dependence of the various neutron 
star properties on the neutron-skin thickness are merely 
due to the fact that different models differ only in the 
density dependence of the symmetry energy. The EOS 
for the symmetric nuclear matter is taken to be the same 
for all the models, since, our goal is to study the diversi¬ 
ties in the properties of the neutron stars arising purely 
due to the differences in the density dependence of the 
symmetry energy within the extended RMF model. In 
Fig. [8] the radii and red shifts for the neutron stars with 
mass 1.4Mq are plotted against the Arnp. The spread in 
i?i ,4 and ZiA for several cases corresponding to the simi¬ 
lar Arnp within the same family is smaller. The values of 
i?i ,4 and Zi .4 obtained for two different families differ no¬ 
ticeably at the smaller values of Arnp. For Arnp = 0.15 
















FIG. 8: (Color online) Plots for the radius R\a for the neutron 
stars with the canonical mass 1.4M0 (upper panel) and the 
corresponding red shift (lower panel) as a function of Arnp 
in the nucleus obtained for the Fp and F2, families of 

models. 


fm the maximum differences in the values of Ria and 
ZiA obtained for the two families are ^ 1.0 km and 0.02, 
respectively. In Fig. [5] We display our results for thresh¬ 
old mass the M^u required for the enhanced cooling of 
neutron stars by means of neutrino emission from the 
nucleons in the direct Urea process [53|. The values of 
Mou are quite sensitive to the neutron-skin thickness. 
The value of M^u for the F 2 p family can vary over the 
range of 0.8 — I.OMq with Acnp decreasing from O.Sfm 
to O.lfm. This variation is little smaller for the case of 
Fp family. The value of M^u for both the families differ 
quite significantly at smaller Arnp. At Acnp =0.15 fm 
the difference between the values of M^u for both the 
families is about 0.6Mq which is quite significant (40%). 

We now consider our results for the tidal polarizability 
parameter A defined as, 

Q,j=-X£,j, ( 10 ) 

where, Qij is the induced quadrupole moment of a star 
in binary due to the static external tidal field of the com¬ 
panion star £ij. The parameter A can be expressed in 
terms of the dimensionless quadrupolar tidal Love num¬ 
ber ^2 as, 

A = (11) 

where, R is the radius of a isolated neutron star, i.e., 
long before merger. The value of ^2 depends on the stel¬ 
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FIG. 9: (Color online) The dependence of the threshold neu¬ 
tron star mass Mdu on the neutron-skin thickness in ^°®Pb 
nucleus. The neutron stars with mass equal to or larger than 
Mdu undergo enhanced cooling through direct Urea process 
for the cooling. 


lar structure and can be calculated by following the pro¬ 
cedure used in Refs. [5§ - l^ . The values of A for the 
neutron stars with masses ^ lAT© are sensitive to the 
behaviour of the symmetry energy at supra-nuclear den¬ 
sities [^. In Fig. [TUI ■we plot the values of A as a function 
of neutron star mass obtained for different parameteriza- 
tions for the Fp and F2p families. The value of Arnp is 
equal to 0.15 fm for all of these cases. The differences 
in the tidal polarizability at low mass neutron star for 
the two different families is very small. But the differ¬ 
ence increases as the mass increases due to different high 
density behaviour of the symmetry energy for different 
families of models. The values of A for the neutron star 
with canonical mass vary over a wide range of 2.7 x 10^® 
to 4.3 X 10®® cm^gs^. The value of A at 1.4M© obtained 
for the Fp family is about 1.5 times larger than that for 
the F 2 p family. The inclusion of S mesons slightly low¬ 
ers the value of tidal polarizability of neutron star with 
mass 1.4Mq. On passing, it may be remarked that the 
differences in A across the two different families of the 
models are larger than the uncertainties expected in its 
measurement by the advanced LIGO-’Vergo detector and 
Einstein Telescope [13, [o^l ■ 

Finally, we have collected in Table IITTI the results for 
the various properties of the neutron stars obtained for 
a few representative cases corresponding to the Fp and 
F 2 p families. The values of neutron-skin thickness for the 
208pb nucleus are also listed. The comparison of the re¬ 
sults obtained for the SETl with those for SET3 or SET2 
with SET4 readily gives a crude estimate about the vari¬ 
ation in the properties of the neutron stars across the 
different families of models at a fixed neutron-skin thick¬ 
ness. Similarly, the idea about the effects of 6 mesons 





9 



FIG. 10: (Color online) Variations in the tidal polarizabil¬ 
ity parameter A with the neutron star mass for the different 
parameterizations of Fp and F 2 p families corresponding to 
neutron-skin thickness Arnp = 0.15 fm in the nucleus. 


TABLE III: Properties of neutron stars and the neutron-skin 
thickness in the ^°®Pb nucleus obtained for SETl -SET4 pa¬ 
rameters. The values of tidal polarizability parameter Ai, 4 , 
listed in the last row, correspond to the neutron star with 
mass 1.4M0. 


Properties 

4 


F2p 

SETl 

SET2 

SETS 

SET4 

Arnp (fm) 

0.22 

0.15 

0.22 

0.15 

pDu (fm)“® 

0.297 

0.282 

0.401 

0.505 

pt(fm)“® 

0.058 

0.069 

0.073 

0.107 

Ft (MeV fm"®) 

0.222 

0.107 

0.474 

0.509 

i?i.4(km) 

13.08 

12.96 

13.00 

12.37 

-Hmax (km) 

11.40 

11.68 

11.29 

11.28 

Mdu{Mq) 

1.01 

1.09 

1.33 

1.69 

Ai, 4 (10^® cm^g s^) 

3.41 

4.33 

2.87 

2.88 


within the same family can be obtained by comparing 
the results for SETl with SET2 or SETS with SET4. It 
may be easily verified from Table Hill and Figs. [8] - [10] 
that the values of core-crust transition density pt,RiA, 
Mou and the tidal polarizability parameter A for both 
the families of the models can differ significantly at a 
fixed value of Ar^p. Thus, instead of the a — p and 
uj — p cross-couplings as included separately in the dif¬ 
ferent families of the models, a linear combinations of 
these cross-couplings in a single model would allow one 
to adjust the properties of the neutron stars over a wide 
range at a fixed value of the neutron skin thickness in a 
heavy nucleus, like, ^'’®Pb. Furthermore, the presence of 
d mesons enable ones to obtain the models with smaller 
value of neutron-skin thickness as can be seen from Fig. 

El 


V. CONCLUSIONS 

We have studied the differences in the various proper¬ 
ties of the neutron stars arising mainly due to the uncer¬ 
tainties in the density content of the nuclear symmetry 
energy in the extended RMF model. With this aim, two 
different families of the extended RMF model, namely, 
Fp and F 2 p are obtained. The Fp family includes a — p 
cross-coupling, while, the F 2 p family includes lo — p cross- 
coupling. Both the families of models include the contri¬ 
butions from the 5 meson in addition to several linear and 
non-linear interaction terms already present in the com¬ 
monly used RMF models. The several parameterizations 
for each of the families of the models are so obtained that 
they yield wide variations for the neutron-skin thickness 
Arnp in the ^°®Pb nucleus without affecting much the 
other bulk properties of the nuclei. The inclusion of 5 
meson produces required splitting in the effective mass 
for the neutrons and protons and also enables us to ob¬ 
tain smaller neutron-skin thickness. 

We compare the various properties of neutron stars 
obtained for the Fp and Fjp families. The properties 
of neutron stars considered are the core-crust transition 
density, radius, red shift, tidal polarizability parameter 
and threshold mass required for the enhanced cooling 
through direct Urea process. Most of these properties of 
the neutron stars at a fixed Arnp are noticeably different 
for two different families of the models. These differences 
are pronounced at smaller values of Arnp which can be 
attributed to the differences in the density dependence 
of the symmetry energy resulting due to different cross¬ 
coupling terms. For Arnp = 0.15fm in the ^°®Pb nucleus, 
consistent with the current experimental data on dipole 
polarizability , the red-shift and the radius of neutron 
stars with mass 1.4Mq differs by about 10% for the two 
families of models. Such differences are quite significant 
(~ 40%) for the tidal polarizability parameter for the 
neutron stars with mass 1.4Mq and the threshold mass 
required for the direct Urea process to occur in the neu¬ 
tron stars. The values of the core-crust transition density 
also differs reasonably across the different families of the 
models. We may thus say that the simultaneous inclusion 
of the a — p and w — p cross-couplings in the extended 
RMF model would enhance its flexibility to accommo¬ 
date the variations in the properties of the neutron stars 
at a given neutron-skin thickness. 
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